Spin-orbit coupling in a graphene bilayer and in graphite 
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The intrinsic spin-orbit interactions in bilayer graphene and in graphite are studied, using a tight 
binding model, and an intraatomic LS coupling. The spin-orbit interactions in bilayer graphene and 
graphite are larger, by about one order of magnitude, than the interactions in single layer graphene, 
due to the mixing of 7r and a bands by interlayer hopping. Their value is in the range 0.1 — IK. The 
spin-orbit coupling opens a gap in bilayer graphene, and it also gives rise to two edge modes. The 
spin-orbit couplings are largest, ~ 1 — 4K, in orthorhombic graphite, which does not have a center 
of inversion. 



INTRODUCTION 

The isolation and control of the number of carriers 
in single and few layer graphene flakes [U [2] has lead to 
a large research activity exploring all aspects of these 
materials [3]. Among others, the application of graphene 
to spintronic devices [3HTU] and to spin qubits [TTHT5] is 
being intensively studied. The understanding of these 
devices requires a knowledge of the electronic spin-orbit 
interaction. In principle, this interaction turns single 
layer graphene into a topological insulator [2], which 
shows a bulk gap and edge states at all boundaries. 
The magnitude of the spin-orbit coupling in single layer 
graphene has been studied [l5Tfl"8] . The calculated cou- 
plings are small, typically below 0.1K. The observed spin 
relaxation [5J[T2] suggests the existence of stronger mech- 
anisms which lead to the precession of the electron spins, 
like impurities or lattice deformations |20f[2"2"] . 

Bilayer graphene is interesting because, among other 
properties, a gap can be induced by electrostatic means, 
leading to new ways for the confinement of electrons |23j. 
The spin-orbit interactions which exist in single layer 
graphene modulate the gap of a graphene bilayer [21]. 
The unit cell of bilayer graphene contains four carbon 
atoms, and there are more possible spin-orbit couplings 
than in single layer graphene. 

We analyze in the following the intrinsic and extrinsic 
spin-orbit couplings in bilayer graphene, using a tight 
binding model, and describing the relativistic effects 
responsible for the spin-orbit interaction by a LS in- 
traatomic coupling. We use the similarities between the 
electronic bands of a graphene bilayer and the bands of 
three dimensional graphite with Bernal stacking to gen- 
eralize the results to the latter. 



THE MODEL 

. We describe the electronic bands of a graphene bi- 
layer using a tight binding model, with four orbitals, the 
2s and the three 2p orbitals, per carbon atom. We con- 
sider hoppings between nearest neighbors in the same 



plane, and nearest neighbors and next nearest neighbors 
between adjacent layers, see [25] . The couplings between 
each pair of atoms is parametrized by four hoppings, 
V ss , V sp , V pp7r and V ppa . The model includes also two in- 
traatomic levels, e s and e p , and the intraatomic spin-orbit 
coupling 

W-so = Ago LjSj (1) 

i 

The parameters used to describe the n bands of 
graphite [261 [27], 7o, 7i>72,73, 74, 75 and A, can be de- 
rived from this set of parameters. We neglect the differ- 
ence between different hoppings between atoms which are 
next nearest neighbors in adjacent layers, which are re- 
sponsible for the difference between the parameters 73 
and 74. We also set the difference in onsite energies 
between the two inequivalent atoms, A to zero. The 
parameters 72 and 75 are related to hoppings between 
next nearest neighbor layers, and they do not play a role 
in the description of the bilayer. The total number of 
parameters is 15, although, without loss of generality, 
we set t p = 0. We do not consider hoppings and spin 
orbit interactions which include d levels, although they 
can contribute to the total magnitude of the spin-orbit 
couplings[I51 12"8], The effects mediated by d orbitals do 
not change the order of magnitude of the couplings in 
single layer graphene, and their contribution to interlayer 
effects should be small. 

The main contribution to the effective spin-orbit at 
the Fermi level due to the interlayer coupling is due to 
the hoppings between p orbitals in next nearest neigh- 
bor atoms in different layers. This interaction gives rise 
to the parameters 73 and 74 in the parametrization of 
the bands in graphite. For simplicity, we will neglect 
couplings between s and p orbitals in neighboring layers. 
The non zero hoppings used in this work are listed in 
Table H 

The hamiltonian can be written as a 32 x 32 matrix 
for each lattice wavevector. We define an effective hamil- 
tonian acting on the 7r, or p z , orbitals, by projecting out 
the rest of the orbitals: 

u e J f =n^+n v<T {u-n aa )- x n^ (2) 
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TABLE I: Non zero tight binding parameters, in eV, used in 
the model. The hoppings are taken from[29, 30'., and the spin- 
orbit coupling from[3T]. Superindices 0,1, and 2 correspond to 
atoms in the same layer, nearest neighbors in different layers, 
and next nearest neighbors in different layers. 




FIG. 1: (Color online). Unit cell of a graphene bilayer. La- 
bels A and B define the two sublattices in each layer, while 
subscripts 1 and 2 define the layers. 



We isolate the effect of the spin-orbit coupling by defin- 
ing: 
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FIG. 2: (Color online). Dependence on energy of the spin- 
orbit couplings, as defined in eq. [5] 



These symmetries are defined by the operators: 
T 



--120 





(4) 



where K. is complex conjugation. 

The possible spin dependent terms which respect these 
symmetries were listed in|32j. in connection with the 
equivalent problem of three dimensional Bernal graphite 
(see below). In the notation described above, they can 
be written as 

~W° = \\o z t z s z + A 2 /i z t 2 s z + A 3 ^ z {<J y s x - T z a x s y ) + 

+ A4<7 Z (t-t-ySx + T z /I x S y 

) (5) 

The first term describes the intrinsic spin-orbit coupling 
in single layer graphene. The other three, which involve 
the matrices are specific to bilayer graphene. The 
term proportional to A3 can be viewed as a Rashba cou- 
pling with opposite signs in the two layers. 



RESULTS 



Note that depends on the energy, uj. 

We analyze at the K and K' points. The two 
matrices have a total of 16 entries, which can be labeled 
by specifying the sublattice, layer, spin, and valley. We 
define operators which modify each of these degrees of 
freedom using the Pauli matrices a, (I, s, and f. The unit 
cell is described in Fig. [T] 

The hamiltonian has inversion and time reversal sym- 
metry, and it is also invariant under rotations by 120°. 



Bilayer graphene 

The energy dependence of the four couplings in eq. [5] 
is shown in Fig. [2j The values of the couplings scale 
linearly with A so . This dependence can be understood 
by treating the next nearest neighbor interlayer coupling 
and the intratomic spin-orbit coupling as a perturbation. 
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FIG. 3: (Color online). Dependence on interlayer gap, E g , of 
the spin-orbit couplings, as defined in eq. [5] 

The spin-orbit coupling splits the spin up and spin down 
states of the a bands in the two layers. The interlayer 
couplings couple the 7r band in one layer to the a band in 
the other layer. Their value is of order 73. The tt states 
are shifted by: 

~ -r^-r oc T&so ( ^) (6) 

where e° is an average value of a level in the a band. 

The model gives for the only intrinsic spin-orbit cou- 
pling in single layer graphene the value 

|Af LG | = 0.0065meV (7) 

This coupling depends quadratically on A so , <5e„.± ~ 

The band dispersion of bilayer graphene at low ener- 
gies, in the absence of spin-orbit couplings is given by 
four Dirac cones, because of trigonal warping effects as- 
sociated with 73 [23]. Hence, we must to consider the 
couplings for wavevectors k slightly away from the K 
and K' points. We have checked that the dependence of 
the couplings A, on momentum, in the range where trig- 
onal warping is relevant, is comparable to the changes 
with energy shown in Fig. [2] 

A gap, E g , between the two layers breaks inversion 
symmetry, and can lead to new couplings. The calcula- 
tions show no new coupling greater than 10 _6 meV for 
gaps in the range — O.leV < E g < O.leV. The depen- 
dence of the couplings on the value of the gap is shown 
in Fig. [3] This calculation considers only the effect in the 
shift of the electrostatic potential between the two layers. 
The existence also of an electric field will mix the p z and s 
orbitals within each atom, leading to a Rashba term sim- 
ilar to the one induced in single layer graphene [TBI fTB] . 

The effect of Ai is to open a gap of opposite sign in 
the two valleys, for each value of s z . The system will be- 
come a topological insulator |141 155] . The number of edge 
states is two, that is, even. The spin Hall conductivity 



FIG. 4: (Color online). Dependence on momentum perpen- 
dicular to the layers in Bernal graphite of the spin-orbit cou- 
plings, as defined in eq. [5] 

is equal to two quantum units of conductance. A per- 
turbation which preserves time reversal invariance can 
hybridize the edge modes and open a gap. Such pertur- 
bation should be of the form t x s v . 

The terms with A3 and A4 describe spin flip hoppings 
which involve a site coupled to the other layer by the 
parameter 71. The amplitude of the wavefunctions at 
these sites is suppressed at low energies [25]. The shifts 
induced by A3 and A4 in the low energy electronic levels 
will be of order A3/71, A4/71. 

Bulk graphite 

. The hamiltonian of bulk graphite with Bernal stack- 
ing can be reduced to a set of bilayer hamiltonians with 
interlayer hoppings which depend on the momentum 
along the direction perpendicular to the layers, k z . We 
neglect in the following the (small) hoppings which de- 
scribe hoppings between next nearest neighbor layers, 72 
and 75 , and the energy shift A between atoms in different 
sublattices. At the K and K' points of the three dimen- 
sional Brillouin Zone (2k z c = 0, where c is the interlayer 
distance) the hamiltonian is that of a single bilayer where 
the value of all interlayer hoppings is doubled. At the H 
and H' points, where 2k z c = 7r, the hamiltonian reduces 
to two decoupled layers, and in the intermediate cases the 
interlayer couplings are multiplied by |2cos(fc z c)|. Carry- 
ing out the calculations described in the previous section, 
k z dependent effective couplings, X i (k z ) 1 can be defined. 
These couplings are shown in Fig. [4] The results for 
bilayer graphene correspond to k z c = 27r/3,47r/3. The 
layers are decoupled for k z c = tt. In this case, the only 
coupling is Ai, which gives the coupling for a single layer, 
given in eq. [7] 

The significant dispersion as function of momentum 
parallel to the layers shown in Fig. [4] implies the exis- 
tence of spin dependent hoppings between layers in differ- 
ent unit cells. This is consistent with the analysis which 



4 



lAil(meV) 




|A 2 |(meV) 




FIG. 5: (Color online), 
the spin-orbit couplings for orthorhombic graphite, as defined 
in eq. [9] The point k x = 0, k y a\/3 = 47r/3 corresponds to the 
K point (a is the distance between carbon atoms in the plane). 



showed that the spin-orbit coupling in a bilayer has a 
contribution from interlayer hopping, see eq. [6] 

The spin-orbit couplings can be larger in bulk graphite 
than in a graphene bilayer. The bands in Bernal graphite 
do not have electron-hole symmetry. The shift in the 
Fermi energy with respect to the Dirac energy is about 
Ep k, 20meV Ai,A3[3I]. Hence, the spin-orbit cou- 
pling is not strong enough to open a gap throughout the 
entire Fermi surface, and graphite will not become an 
insulator. 

A similar analysis applies to orthorhombic graphite, 
which is characterized by the stacking sequence 
ABC ABC ■ ■ ■ [35]. The electronic structure of this al- 
lotrope at low energies differs markedly from Bernal 
graphite[36l [37], and it can be a model for stacking 
defects [36H38] , If hoppings beyond nearest neighbor lay- 
ers are neglected, the hamiltonian can be reduced to an 
effective one layer hamiltonian where all sites are equiv- 
alent. The effective hamiltonian which describes the K 
and K' valleys contains eight entries, which can be de- 
scribed using the matrices <Xi } s-j, and Ti. Orthorhombic 
graphene is not invariant under inversion, and a Rashba 
like spin-orbit coupling is allowed. The spin-orbit cou- 
pling takes the form: 



u 



SO 

ortho 



\ortho _ _ 
A x O z S z T z 



\ ortho { _ 

A 2 {a y 



T z (J x Sy 



(8) 



As in the case of Bernal stacking, the couplings have a 
significant dependence on the momentum perpendicular 
to the layers, k z , and interlayer hopping terms are in- 
duced. For to = 0, k = and k z = 0, we find: 



K 

AS 



ortho 



rtho 



0.134meV 
0.275meV 



(9) 



In orthorhombic graphite the Fermi level is away from 
the K and K' points, in the vicinity of a circle defined 
by |k| = 7! /v F |3j)l [37] . The variation of the couplings as 
function of wavevector is shown in Fig. [5] 



CONCLUSIONS 

We have studied the intrinsic spin-orbit interactions in 
a graphene bilayer and in graphite. We assume that the 



origin of the couplings is the intraatomic LS interaction, 
and we use a tight binding model which includes the 2s 
and 2p atomic orbit als. 

The intrinsic spin-orbit couplings in a graphene bilayer 
and in graphite are about one order of magnitude larger 
than in single layer graphene, due to mixing between 
the 7r and a bands by interlayer hoppings. Still, these 
couplings are typically of order 0.01 — O.lmeV, that is, 
0.1 - IK. 

Bilayer graphene becomes an insulator with an even 
number of edge states. These states can be mixed by per- 
turbations which do not break time reversal symmetry. 
These perturbations can only arise from local impurities 
with strong spin-orbit coupling, as a spin flip process and 
intervalley scattering are required. 

The interplay of spin-orbit coupling and interlayer hop- 
ping leads to spin dependent hopping terms. The spin- 
orbit interactions are largest in orthorhombic graphite, 
which does not have inversion symmetry. 
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